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C^ : Abstract 

It is argued, that adjusting strong potentials directly to observed hadronic atom level shifts niay 
lead to significantly different scattering lengths, than those, predicted by the Deser formula l|. 
On the example of the Is level shift of kaonic hydrogen it is demonstrated, that the usually 
adopted Deser jyalues deduced from the two recent measurements in KEK [^ and by the DEAR 



CN ■ Collaboration [3] aD(KEK) = 0.78 - 0.49z fm and aD(DEAR) = 0.47 - 0.3i fm should be replaced 

> 
^^ \ by as(KEK) ~ 0.85 - 0.62i fm and as(DEAR) ~ 0.49 - 0.35z fm, correspondingly. 
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I. INTRODUCTION 

One of the most relevant sources of information on hadron-nucleus interaction is the 
measurement of hadronic atom level shifts. The usual picture of a hadronic atom is a 
particle moving in the combined Coulomb and nuclear potential of the nucleus, (see e.g. |^). 
The measured level shifts are used to extract information about the nuclear part of the 
potential. 

An important subfield is the study of hadronic hydrogen, which gives information on 
the primary hadron nucleon interaction. According to the usual philosophy of hadronic 
atoms, this interaction is imagined as a complex potential, the properties of which are to be 
deduced (at least partly) from the measured level shifts. On the other hand, the elementary 
hadron nucleon interaction can be approached also from the field-theoretical side, trying 
to derive it from effective Lagrangians and then to relate it to the level shifts p]. While 
the field theoretical approach can be considered as more fundamental in this case, the basic 
merit of the potential picture leading to a phenomenological potential lies in its applicability 
in dynamical description of more complicated N > 2 systems in the quantum mechanical 
framework. Proper field theoretical calculations of such systems are still beyond the real 
possibilities. 

In both cases the level shifts and the nuclear interactions are related via the scattering 
length of the nuclear interaction. The use of scattering length seems to be natural in 
the field-theoretical case since it is the zero energy scattering amplitude and thus diagram 
technique can be applied for its calculation. On the other hand, in the potential approach 
this intermediate step seems to be superfluous, with present day computational facilities the 
potentials can be directly and easily related to the level shifts. 

The main purpose of the paper is to show, that the approximate deduction of the potential 
parameters via the widely used Deser- formula |1], connecting the level shift and the scat- 
tering length can lead to substantial differences, compared to the direct approach (this has 
been known before). We argue, that the effort of exact solution of the Schrodinger-equation, 
yielding accurate potential parameters is not greater, than that of applying different correc- 
tions to the Deser-formula. Also, a new approximate relation between the level shifts and 
the potentials is derived, which practically gives the exact results. 

As a demonstration of the above ideas we have made calculations for the the kaonic 



hydrogen case since recently the KN interaction has attracted some renewed interest due 
to the possible existence of bound ^-nuclear states. 

II. FORMULATION OF THE PROBLEM 

Let us consider the Is level of kaonic hydrogen. The (model) Hamiltonian reads 

/7 = -^Ar + K + K, (1) 

where the strong interaction is represented by a local central potential 

Vs = V,Vs{r/h) (2) 

with two parameters: its strength Vq and range h. The attractive Coulomb potential is 

Vc = — (3) 

r 

and yU is the K~p reduced mass: 

/x = 323.478 MeV/cl (4) 

The Schrodinger equation {H — E)'^ = Q can be transformed into the radial equation 

fd' .2 , .. ,. ..^ , 2/ro 

where 



^^2 g^-At;.(r/6) + ^)^(r) = 0, (5) 



and tq is the Bohr radius 

ro = ^ = 83.594 fm. 

The "atomic" energy unit in this case is 

4 

60 = ^ = 17.226keV. 

In connection with Eq. ([5]) we are interested in the change of the eigenvalue g^ compared 
to the Is eigenvalue Qq of the pure Coulomb equation 

rf2 2 , 2/^0 



dr2 -ro + ^] V{r) = 0. (6) 



The solution of Eq. (^ is, of course, known: 

1 r f r 

qo = —, ^ exp 

ro ro V '^o 

Let us recall some important features of Eq. ([5]): 



— the strong interaction acts on the fm scale and is vanishingly small beyond a certain 

distance i?: 6 ~ fm, Vs{r/b) ~ for r > i?; 

— due to absorbtion to other channels, A is complex, so q is complex, too, but Re{q) > 0, 

so il){r) —>■ for r —>■ 0; 

— both functions iplr) and (p{r) are on the tq scale (~ 100 fm), while the presence of Vs 

modifies ip{r) compared to (p{r) only for r < R; thus g^ — q^ is small, but not due to 
the smallness of Vg itself, so perturbative treatment is not justified; 

— outside the range of the strong potential (r > R) ipijr) goes over into the asymptotically 

vanishing solution of the Coulomb equation 

'^i.'r) ~ KuM^'r) = rexp{-qr)U i 1 , 2; 2gr 

with U being the confluent hypergeometric function of the second kind 6(]. 

— within the range of Vs (r < R) the relation g^ <^ Xvs{r/b) usually holds since the 

Coulombic eigenvalue is in the keV range, while the nuclear potential is of MeV order; 
this feature allows to approximate ^(r) in this range by the zero-energy solution ^°(r). 

Multiplying Eq. ([5]) by ip{r) and Eq. (j6]) by 4'{r), subtracting and integrating, we obtain 

in the usual way: 

2 _ 2 ^ J^v{r)Xvs{r/b)ij{r)dr 
J^ (p{r) ijj{r) dr 

where we used the fact, that both ilj{r) and !f{r) vanish for r ^ oo and r — > 0. 

The expression ((Tj) is exact and independent of the normalization of ilj{r) and f{r). The 

basic question is how to relate 

totheproperfesofV,. 

The traditional answer to this question was given by Deser p|, back in 1954. According 
to him 



,2 „2 



a.. 



q — qo ^ —4: ^ or 



r, 
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«D(g) = -j(g'-go)~«s, (8) 



where a^ is the scattering length of the strong potential Vs defined as ^ 



' 1 + /J" A i;,(r/6) ^'^^{r) dr €'(^) ' 

with ip^{r) being the regular [V's(O) = 0!"^° (0) = 1] solution of the zero-energy equation 



with K alone: 



s 



2 



l^-Xv,{r/b)jrs{r)=0. (10) 

In other words, the pure strong scattering length is approximated by the Deser scattering 
length aoiq), derived from the measured AE. 

It is not straightforward to relate ([8]) and (Q to the exact expression ([7]) or to point out 
clearly the approximations leading from (I7j) to ([H]), together with criteria for their applica- 
bility (apart from the obvious R <^ vq). 

Since 1954 considerable effort has been devoted to re-derivation of Deser's result, to 
considering its possible improvements or corrections to it, to studying its special cases, e.g. a 
strong potential with an almost or weakly bound state, when the scattering length becomes 



large [7|]-ll|. The most relevant improvement is the taking into account the Coulomb 



distortion of the zero-energy wave function ip^, leading to approximation 

aoiq) ~ ttsc (11) 

instead of ([8]), where age is the Coulomb modified scattering length: 



1 + /o e(r)A Vsir/b) ^Pl{r) dr WiQ, ^oj|,= 



■R 



Here <I>(r) and 0(r) are the suitable q ^ limits of the Coulomb scattering functions F 
and G (see |l2|), satisfying 

£-^){r:Ho (13) 

and ip^ci^) is the regular zero-energy solution of 

^-At^.(r/6) + ^)^°,(r) = 0, (14) 



^ This is the usual textbook definition of as corresponding to negative scattering length for weak attractive 
potential (A < 0). By some reason in meson- nuclear physics the opposite sign is used. 



while W denotes the wronskian of the two functions. 

Now, the usual way, how ([8]) and flTTj) are used to relate a strong potential K to the 
measured energy shift is to approximate a^ (or Qsc) by aoiq) and then to design a potential 
Vg having this a^ (or age)'- 

asc ~ aoiq) 

But due to the approximations ([H]) or flTTl) the A£" values, calculated with the potentials 
obtained in this way do not reproduce the measured AE-s : AE' ^ AE, so the usual 
experimental claim, that measuring AE is equivalent to measuring a^ is not valid. 

However, at the present calculational level (2007 ^ 1954) to find a potential, giving a 
prescribed a^ (or asc) is not easier, than to find a potential giving exactly the measured AE. 
The eigenvalue equation i^ is solved by numerical integration in the internal region (r < R) 
and by matching the logarithmic derivatives of ipinir) and the external function F^^^{q,r) 
at r = R: 

For a given range parameter b the root in q of Eq. flTSj) for fixed A gives the eigenvalue, 
while the root in A for fixed q yields the potential strength corresponding to a prescribed 
eigenvalue. The numerical solution of Eq. flT^ is straightforward in both cases. 

For demonstration of the above idea strong potentials were derived from the three con- 
ditions 

AE(calc) = AE(exp) (a) 

tts = aD(gexp) (b) 

asc = aD(gexp) , (c) 

their level shifts and strong scattering lengths were compared for four type of commonly used 
potential shapes: exponential, gaussian, square well, and Yamagouchi (non-local, separable) 
and different ranges b. The experimental AE{ex'p) was taken from the KEK experiment [2]: 

AE(KEK) = -323 + 203.5z eV with 

aoiKEK) = 0.78 - 0.49i fm . (16) 



III. RESULTS AND CONCLUSIONS 

The main message of this work is to emphasize, that no approximations are needed to 
connect the measured A£^(exp) with the properties of the strong interaction (model) poten- 
tial Vg. However, the desire to find a relation between these quantities without solving the 
eigenvalue equation, which is superior to the previously used ones, motivated the derivation 
of another approximate formula. 

The normalization integral in the denominator of Eq. ([7]) can be approximated as 

oo ^oo 

(f{r)-ip{r)dr^ / ip{r)tp{r) dr = (17) 

Jr 



poo 

Jr 



due to the smallness of the nuclear region < r < R compared to the whole range of func- 
tions f{r) and ip{r). Using the equations satisfied by f{r) and F^^^{r) and their vanishing 
for r ^ oo Eq. (IT71) can be rewritten as 

■fir)FUl.r)dr^^^^il4^. (18) 

R q^ -Qo 

Substituting (fTSjl into (TT71) and ([7j) and, as before, approximating ip{r) within the nuclear 
range by ^scl'") ^^ get finally 

W{F:^,,^)\r=R =- f v{r)\vs{r/h)^l{r)dT , (19) 

which is of the form 

w{q) = /(A) (d) 

and again can be solved either for q or for A. The results for the strong potentials obtained 
from condition (d) are also shown in Table 1, summarizing our results. 
From these results we can make the following conclusions: 

(i) The quality of approximations 

It can be seen, that the most commonly used approximation (6) gives rather poor 
results, in the considered case the error in AE' is of the order of 10 — 15%. The 
opposite is also true: the a^-s corresponding to potentials exactly reproducing the 
measured AE differ from the Deser value an approximately by the same amount. The 
improved approximation (c) yields better strong potentials, the error both in AE' and 



Qs amounts to a few %. The best results seem to be given by approximation ((i), which 
practically reproduces the exact ^E-s and a^-s of the exact calculation (a). 

(ii) Model independence 

One of the most attractive features of the Deser formula is its model independence: the 
relation between the energy shift and the scattering length is independent of the form 
of the potential. This feature is confirmed by the calculations with different potential 
shapes and quite different ranges. The relation between AE and a^ is almost potential 
independent, however, only within one approximation. 

(iii) Finally, to the question ,, which feature of a given potential instead of as or asc deter- 
mines AE ? ", our answer is: 

rR 

/(A) = - / i^{r)Xvs{r/h) ^°,(r) dr 
Jo 

through the relation (d). 

(iv) The widely used ,, measured" values of as(KEK) for the kaonic hydrogen should be 
changed from a£)(gexp) (IT6l) to 

a,(KEK) ^ 0.85 - 0.62z fm 

since this is the value, which corresponds to potentials exactly reproducing the mea- 
sured level shift. 

The same calculation was also performed for the case of the DEAR Is level shift [3|, 
and while the overall results are qualitatively the same as before, the correct scattering 
length corresponding to the DEAR shift is 

a, (DEAR) ~ 0.49 - 0.35i fm 

instead of the widely adopted a£)(DEAR) = 0.47 — 0.30i fm. 
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TABLE I: Calculated kaonic hydrogen Is level shifts AE and scattering lengths Os of dif- 
ferent strong interaction potentials, derived from the measured KEK level shift A£^(exp) = 





-323 + 203.5i eV using 


the conditions (a) — {d] 


1 (see the text). 






Gauss 


Exponential 


Square well 


Yamaguchi 


b, fm 


as, fm AE, eV 


Os, fm AE, eV 


as, fm AE, eV 


as, fm AE, eV 


Vs from AE(calc)=AE(exp) 



1.0 0.85-0.61i 0.88-0.61i 0.84 - 0.61i 0.86 - 0.61i 

0.5 0.84 - 0.62i -323 + 204i 0.85 - 0.6U -323 + 204i 0.84 - 0.62i -323 + 204i 0.84 - OMi -323 + 204i 
0.25 0.84-0.63i 0.84 - 0.62i 0.84 - 0.64i 0.84 - 0.62i 

Vs from as = ao 

1.0 -298 + 168i -289 + 167i -299 + 167i -294 + 168i 

0.5 0.78 - 0A9i -299 + 165i 0.78 - 0A9i -296 + 168i 0.78 - 0.49i -299 + 164i 0.78 - 0A9i -298 + 167i 
0.25 -299 + 162i -298 + 166i -298 + Wli -300 + 165i 

Vs from Use = ttD 

1.0 0.83 - 0.57i -317 + 191i 0.86 - 0.57i -317 + 19li 0.83 - 0.57i -317 + Wli 0.84 - 0.58i -317 + 19H 
0.5 0.83 - 0.58i -317 + Wli 0.84 - 0.57i -317 + 191i 0.83 - 0.58i -317 + Wli 0.83 - 0.57i -317 + 191i 
0.25 0.83 - 0.59i -317 + 191i 0.83 - 0.58i -317 + Wli 0.83 - 0.60i -317 + Wli 0.83 - 0.56i -317 + 19H 



T4 from w{q) = /(A) 



1.0 0.85 - OMi -323 + 204i 0.86 - 0.61i -319 + 202i 0.84 - 0.61i -323 + 20Ai 0.86 - OMi -323 + 204i 
0.5 0.84 - 0.62i -323 + 20Ai 0.86 - OMi -319 + 202i 0.84 - 0.62i -322 + 204i 0.84 - 0.61i -322 + 203i 
0.25 0.84 - OMi -323 + 20Ai 0.85 - 0.61i -322 + 203i 0.84 - 0.63i -323 + 20Ai 0.84 - 0.62i -323 + 203i 



